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The thermal multihadron production observed in different high energy colhsions poses two basic 
problems: (1) why do even elementary collisions with comparatively few secondaries (e"'"e~ annihi- 
lation) show thermal behaviour, and (2) why is there in such interactions a suppression of strange 
particle production? We show that the recently proposed mechanism of thermal hadron production 
through Hawking-Unruh radiation can naturally account for both. The event horizon of colour con- 
finement leads to thermal behaviour, but the resulting temperature depends on the strange quark 
content of the produced hadrons, causing a deviation from full equilibrium and hence a suppres- 
sion of strange particle production. We apply the resulting formalism to multihadron production 
in e^e^ annihilation over a wide energy range and make a comprehensive analysis of the data in 
the conventional statistical hadronization model and the modified Hawking-Unruh formulation. We 
show that this formulation provides a very good description of the measured hadronic abundances, 
fully determined in terms of the string tension and the bare strange quark mass; it contains no 
adjustable parameters. 



I. INTRODUCTION 



Hadron production in high energy collisions shows remarkably universal thermal features. In e^e~ 
annihilation 3, 0) @]) PP^ PP \M ^^'^ more general hh interactions [3], as well as in the collisions of 
heavy nuclei 5|, over an energy range from around 10 GeV up to the TeV range, the relative abundances 
of the produced hadrons appear to be those of an ideal hadronic resonance gas at a quite universal 
temperature Th ~ 160 — 170 MeV, as illustrated in Fig. [T] The transverse momentum spectra of 
the hadrons produced in hadronic collisions at intermediate energies are also in good agreement with the 
predictions of a statistical model based on the same temperature [3|. 

There is, however, one important non-equilibrium effect observed: the production of strange hadrons 
in elementary collisions is suppressed relative to an overall equilibrium. This is usually taken into account 
phenomenologically by introducing an overall strangeness suppression factor 7s < 1 0], which reduces 
the predicted abundances by 7^,7^ and 7^ for hadrons containing one, two or three strange quarks (or 
antiquarks), respectively. In high energy heavy ion collisions, strangeness suppression becomes less and 
disappears at high energies [8(. 

When comparing the temperatures of different collision channels in Fig. [U it should be stressed that in 
elementary collisions, in contrast to heavy ion collisions, the conservation of charge, baryon number and 
strangeness is enforced exactly (canonical ensemble), so that here no chemical potential is introduced. The 
effect of exact charge conservation is important in elementary collisions because of the low multiplicities 
involved, while it is generally negligible in large multiplicity high energy nuclear reactions. The fitted 
values of 7s lie around 0.5 to 0.7 in elementary collisions; the heavy ion values appear to vary as a function 
of energy, but it tend to approach unity \§\ . 

The origin of the observed thermal behaviour has been an enigma for many years and there is a 
still ongoing debate about the interpretation of these results 0. While the common belief is that in 
high energy heavy ion collisions multiple parton scattering could lead to kinetic thermalization through 
multiple scattering, e~^e~ or elementary hadron interactions do not readily allow such a description. The 
universality of the observed temperatures, on the other hand, suggests a common origin for all high energy 
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FIG. 1: Hadronization temperatures (from ref. 



collisions, and it was recently proposed [T3| that thermal hadron production is the QCD counterpart of 
Hawking-Unruh radiation [lll.ll2l|. emitted at the event horizon due to colour confinement. A well-known 
instance of this phenomenon is the Schwinger mechanism [l3L Il4l | of pair production in a constant electric 
field £. The probability of spontaneously producing an electron-positron pair is in this case given by 

P{m.£) ^ exp{-7rm^/e£}, (1) 

where m is the electron mass and e its charge. Since 

2e£ 



m 



(2) 



is just the acceleration of an electron (of reduced mass to/2) in the field £, we find that the pair production 
probability has the thermal form 

P(TO,f) -cxp{-TO/rc/}, (3) 

where 

In TTTO 



is the corresponding temperature. In other words, it is just that obtained by Unruh IJ] for the radiation 



emitted when a mass suffers constant acceleration and hence encounters an event horizon. 

In the case of QCD and approximately massless quarks, the resulting hadronization temperature is 
determined by the string tension cr, with T ~ y^a/2Tr. The aim of the present work is to show that 
strangeness suppression occurs naturally in this framework, without requiring a specific suppression 
factor. The crucial role here is played by the non-negligible strange quark mass, which modifies the 
emission temperature for such quarks. 

In this work, we will focus on e^e~ collisions, which is the simplest case for this model to be tested. 
In the next section, we will briefly review the usual formulation of the statistical hadronization model, 
since this will form the basis also for the subsequent description in terms of Hawking-Unruh radiation 
model, to be formulated in Sect. 3. There we shall in particular derive the dependence of the radiation 
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temperature on the mass of the produced quark and show how this affects the hadronization in e^e~ 
annihilation. In Sect. 4, we then present an updated comprehensive analysis of all available data from 
y/s — 14 GeV to ^/s — 189 GeV and compare the results of the Hawking-Unruh formulation to the 
conventional statistical description. 

II. STATISTICAL HADRONIZATION IN e+e" COLLISIONS 

In this section we will briefly review the essentials of the statistical hadronization model and its appli- 
cation to e"^e~ collisions. For a detailed description see ref. p^ . 

The statistical hadronization model assumes that hadronization in high energy collisions is a universal 
process proceeding through the formation of multiple colourless massive clusters (or fireballs) of finite 
spacial extension. These clusters are taken to decay into hadrons according to a purely statistical law: 
every multi-hadron state of the cluster phase space defined by its mass, volume and charges is equally 
probable. The mass distribution and the distribution of charges (electric, baryonic and strange) among 
the clusters and their (fluctuating) number are determined in the prior dynamical stage of the process. 
Once these distributions are known, each cluster can be hadronized on the basis of statistical equilibrium, 
leading to the calculation of averages in the microcanonical ensemble, enforcing the exact conservation of 
energy and charges of each cluster. 

Hence, in principle, one would need the mentioned dynamical information in order to make definite 
quantitative predictions to be compared with data. Nevertheless, for Lorentz-invariant quantities such 
as multiplicities, one can introduce a simplifying assumption and thereby obtain a simple analytical 
expression in terms of a temperature. The key point is to assume that the distribution of masses and 
charges among clusters is again purely statistical .3;], so that, as far as the calculation of multiplicities is 
concerned, the set of clusters becomes equivalent, on average, to a large cluster (equivalent global cluster) 
whose volume is the sum of proper cluster volumes and whose charge is the sum of cluster charges 
(and thus the conserved charge of the initial colliding system). In such a global averaging process, the 
equivalent cluster generally turns out to be large enough in mass and volume so that the canonical 
ensemble becomes a good approximation of the more fundamental microcanonical ensemble fl6|; in other 
words, a temperature can be introduced which replaces the a priori more fundamental description in 
terms of an energy density. 

It should be stressed that in such an analysis of multiplicities, temperature has essentially a global 
meaning and not local as in hydrodynamical models. The only local meaningful quantity is cluster's 
energy density, and even though the globally fitted temperature value is closely related to it, this does 
not mean that single physical clusters can be described in terms of a temperature, unless they are 
sufficiently large (about 10 GeV in mass, see jT6|). 

In the statistical hadronization model supplemented by global cluster averaging, the primary multi- 
plicity of each hadron species j is given by [3| 

^"■'■^ - ^^^^ ~ \~) ZiQ) 

where V is the (mean) volume and T the temperature of the equivalent global cluster. Here ^(Q) is the 
canonical partition function depending on the initial abelian charges Q — [Q, N, S,C, B), i.e., electric 
charge, baryon number, strangeness, charm and beauty, respectively. We denote by rrij and Jj the mass 
and the spin of the hadron j, and — (Qj, Nj, Sj,C'j, Bj) its corresponding charges; 7s is the extra 
phenomenological factor implementing a suppression of hadrons containing strange valence quarks 
(see Sect. 1). In the sum ([S]), the upper sign apphes to bosons and the lower sign to fermions. For 
temperature values of 160 MeV or higher, Boltzmann statistics corresponding to the term n ~ 1 only in 
the series ([SJ is a very good approximation for all hadrons (within 1.5%) but pions. For resonances, the 
formula ([5]) is folded with a relativistic Breit-Wigner distribution of the mass mj . 
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The canonical partition function can be expressed as a multi-dimensional integral 



X exp 



^^(2J, + 1) J d'p log(l±7f-e-V^^/^-"i.-*)±i 



V 



(6) 



where N is the number of conserved abelian charges. Unlike the grand-canonical case, the logarithm 
of the canonical partition function does not scale linearly with the volume. Therefore, the chemical 
factors Z(Q — nq^j)/Z{Q) turn out to be less than unity at finite volume (canonical suppression) and 
only asymptotically reach their grand-canonical limit of unity, for an initially completely neutral system. 
Indeed, they play a major role in determining particle yields in e'''e~ collisions. 

For all energies considered here, the production of heavy flavoured hadrons is negligible in e~^e~ —^ 
qq events, where g is a light quark {u,d,s). As a result, in formula ([5|) the charm and bottom charge 
can be neglected, so that the abelian charge vector reduces to a three-component form Q = (Q, N, S) 
and can be calculated with a numerical integration [J, On the other hand, e~^e~ qq events, 

where q is a heavy quark, always result in the production of two open heavy flavoured hadrons, arising 
from the primary heavy quark-antiquark pair, which subsequently decay into light-flavoured hadrons. In 
the statistical model, this constraint is readily implemented, requiring the number of heavy quarks plus 
antiquarks in the global description to be two; because of the high charm/bottom mass compared to the 
typical temperature value of 160 MeV, the probability of producing extra heavy qq pairs is absolutely 
negligible. The multiplicities of charm (bottom) hadrons in e^e~ — *■ cc {e'^e~ bb) events become in 
the Boltzmann approximation [3, 0| 

where 

^^^^i^Jj + ^)^'TK2[-^), (8) 

and C denotes the canonical partition function as in ([5]), involving only light-flavoured particles, with 
= 3. The indices j, k label charm (bottom) hadrons, while the index i their antiparticles. 

The statistical treatment of heavy quark formation and hadronization as outlined here effectively means 
heavy quarks occur only in e~^e~ — s- qq interactions, leading to one open charm (bottom) hadron and one 
corresponding antihadron per event. The relative rates of the different possible open charm (bottom) 
states, however, are determined by their phase space weights. 



III. STRING BREAKING AND EVENT HORIZONS 



We first outline the thermal hadron production process through Hawking-Unruh radiation for the 
specific case of e'^e~ annihilation (see Fig. The separating primary qq pair excites a further pair 
qiqi from the vacuum, and this pair is in turn pulled apart by the primary constituents. In the process, 
the qi shields the q from its original partner g, with a new qqi string formed. When it is stretched to 
reach the pair production threshold, a further pair is formed, and so on (isl. [l9j. Such a pair production 
mechanism is a special case of Hawking-Unruh radiation ^14| , emitted as hadron qiq2 when the quark 
qi tunnels through its event horizon to become q2- The corresponding hadron radiation temperature is 
given by the Unruh form Th = a/27r, where a is the acceleration suffered by the quark qi due to the 
force of the string attaching it to the primary quark Q. This is equivalent to that suffered by quark q2 
due to the effective force of the primary antiquark Q. Hence we have 
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FIG. 2: String breaking through qq pair production 

where Wq = \JiTi^ + k"^ is the effective mass of the produced quark, with ruq for the bare quark mass and 

kq the quark momentum inside the hadronic system qiqi or g292- Since the string breaks when it 
reaches a separation distance 

Xq^'^^ml + i'Ka/2), (10) 

the uncertainty relation gives us with fcg ~ 1 /xq 

Wq = ,Jml + [o-2/(4m2 + 27ra)] (11) 

for the effective mass of the quark. The resulting quark-mass dependent Unruh temperature is thus given 
by 

T{qq) ^ (12) 

ZlTWq 

Note that here it is assumed that the quark masses for qi and q2 are equal. For niq ~ 0, eq. (fT2| reduces 
to 

r(oo).y^, (13) 

as obtained in fl^. 

If the produced hadron qiq2 consists of quarks of different masses, the resulting temperature has to be 
calculated as an average of the different accelerations involved. For one massless quark (niq ~ 0) and one 
of strange quark mass rus , the average acceleration becomes 

woao + Wsas 2a 
aos = = ■ (14) 

Wo + Ws Wq + Ws 

From this the Unruh temperature of a strange meson is given by 

T{Os) ~ — -, (15) 



with wq ~ ^/l/Qjra. Similarly, we obtain 



T{ss)^- , (16) 



for the temperature of a meson consisting of a strange quark-antiquark pair ((/)). With a ~ 0.2 GeV^, 
eq. dini) gives Tq ~ 0.178 GeV. A strange quark mass of 0.1 GeV reduces this to T(Os) ~ 0.167 GeV and 
T(ss) ~ 157 MeV, i.e., by about 6 % and 12 %, respectively 

The scheme is readily generalized to baryons. The production pattern is illustrated in Fig. [3] and leads 
to an average of the accelerations of the quarks involved. We thus have 

T(OOO) = r(o) ~ — ^ (17) 
2ttwo 



for nucleons, 



for A and E production, 



for ^ production, and 



T(OOs) ^ ^\ (18) 
1-k{Iwq + Ws) 



TlQss) ~ ^ , ^ r (19) 



nsss) = nss) ^ ^ (20) 

for that of fJ's. We thus obtain a resonance gas picture with five different hadronization temperatures, 
as specified by the strangeness content of the hadron in question: T(00) = r(OOO), T{Qs), T{ss) — 
T{sss), T(OOs) and T{Oss). 




FIG. 3: Nucleon (left) and hyperon (right) production in e^e annihilation 



It is important to stress that, in this picture, the primary quarks produced directly in the annihilation 
are not related to Hawking-Unruh radiation, nor are the light quarks with which they eventually combine 
to hadronize. Therefore, leading hadrons (those containing primary quarks) are essentially different 
and should be treated separately from hadrons containing only newly produced quarks. In practice, 
while in the conventional statistical model the same hadronization temperature governs the relative 
probabilities of emitting different species of leading hadrons, in the Hawking-Unruh formulation, we do 
not have any specific prescription for this. The problem of calculating leading hadron yields is relevant 
in e+e^ annihilations because, in contrast to hadronic collisions, the primary heavy quarks c and b 
have large branching ratios and significantly contribute to the production of light-flavoured hadrons 
through the decay chain, especially in the strange sector. Lacking a definite recipe, we chose to calculate 
the relative heavy flavoured hadron yields by using the same temperatures as for light-flavoured ones, 
quoted above, keeping one weight wq fixed and using wq or Ws according to whether the heavy quark 
hadronization occurs through combination with either u, d or with s, respectively. It should be noted out 
that this is not the only option and that different choices may lead to different results. 

The different species-dependent temperatures are to be inserted into the formulae ([5]) and ([7]) of the 
previous section, in order to determine the primary hadron multiplicities. We note at this point a subtle 
conceptual difference between the conventional statistical approach and the Hawking-Unruh formulation. 
The usual statistical description employs, as noted above, a global cluster averaging, with each cluster 
statistically composed. In the Hawking-Unruh scheme, the radiation in each step is a hadron formed 
from the emitted qq pair, not some thermal cluster. Since the hadron can, however, be a highly excited 
resonance, the two descriptions become equivalent in a Hagedorn-type picture proposing resonances made 
up of resonances in a self-similar pattern. 
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The multiplicities obtained in the Hawking-Unruh scheme are, as emphasized, fully determined by the 
two basic parameters of the formulation, the string tension a and the strange quark mass to^. Apart 
from possible variations of these quantities, our description is thus parameter-free. As illustration, we 
show in tableUthe temperatures obtained for u = 0.2 GeV^ and three different strange quark masses. It 
is seen that in all cases, the temperature for a hadron carrying non-zero strangeness is lower than that 
of non-strange hadrons; as we shall show, this leads to an overall strangeness suppression. 
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m.s = 0.075 


= 0.100 


ras = 0.125 


r(oo) 


0.178 


0.178 


0.178 


T(Os) 


0.172 


0.167 


0.162 


T{ss) 


0.166 


0.157 


0.148 


r(ooo) 


0.178 


0.178 


0.178 


T'(OOs) 


0.174 


0.171 


0.167 


T(Oss) 


0.170 


0.164 


0.157 


T(sss) 


0.166 


0.157 


0.148 



TABLE I: Hadronization temperatures for hadrons of different strangeness content, for rus = 0.075,0.100,0.125 
Gev and a = 0.2 GeV^ 



Our picture implies that the produced hadrons are emitted slightly "out of equilibrium" , in the sense 
that the emission temperatures are not identical. As long as there is no final state interference between 
the produced quarks or hadrons, we expect to observe this difference and hence a modification of the 
production of strange hadrons, in comparison to the corresponding full equilibrium values. Once such 
interference becomes likely, as in high energy heavy ion collisions, equilibrium can be at least partially 
restored, weakening the strangeness suppression. The extension of our approach to heavy ion collisions 
will be dealt with in a subsequent paper. 

IV. ANALYSIS OF HADRON MULTIPLICITIES 

Multihadron production in e'^e" annihilation has been studied at PEP, PETRA and LEP over an 
energy range from 14 to 189 GeV, and the multiplicities of a large number of different species have been 
measured. The relevant data and their references are compiled in the Appendix. 

In order to compare the models with the experimental data, we have at each energy made a fit to 
available measured multiplicities of light-flavoured hadrons, both in the traditional statistical model and 
in the Hawking-Unruh formulation. The traditional model has three free parameters to be determined, 
namely the temperature T, the global volume V (see the discussion in the Sect. 2), and the strangeness 
undersaturation parameter 7^. In the Hawking-Unruh model, V is kept, but the string tension a and the 
strange mass rus replace T and 75 as fit parameters. 

In the fit, the theoretical multiplicity of a given species, to be compared to the data, is calculated as 
the sum of primary multiplicity given by Q and the contribution from the decay of unstable heavier 
hadrons, 

(n,) = {n,r'" + ^ Br(A: ^ j){n,), (21) 

k 

where the branc hing ratios are the measured values as listed in the latest compilation of the Review of 
Particle Physics [2(J|. For the decays of heavy flavoured hadrons with unknown branching fractions, we 
have used the predictions of the PYTHIA [2l| program. The hadrons considered unstable in e+e~ ex- 
periments are all species except tt, K^, K^, p, n, and we have followed this convention in the theoretical 
calculation to meet the definition of measured multiplicities. The hadrons and resonances contributing 
to the sum in Eq. (I^T|) consist here of all known states [lO] up to a mass of 1.8 GeV. 

A specific fraction of e^e~ annihilations occurs into heavy c and b quarks. In this case, the multiplicities 
of light-flavoured hadrons are affected by the presence of the heavy quarks, both at primary level as the 
canonical partition function changes (see discussion in the previous section) and at final level because of 
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the heavy flavoured hadron decays. This is taken into account in our calculations, and the production 
rate of the j*'' hadron is given by 

5 

(n,),. (22) 

where j = 1, 5 accounts for u, d, s, c, b quarks and the index qi specifies the corresponding multiplicity 
in e^e~ — > qiiji annihilation. The values of 

R. ^ ;f T f , (23) 
^(c+e hadrons) 

are the corresponding branching fractions, obtained at each centre-of-mass energy from measurements 
and electroweak calculations. We have here taken the values calculated in the PYTHIA programme [2l[, 
which are quoted in table [iTl for each centre-of-mass energy. 
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0.46 


0.09 
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0.36 


0.09 


91.25 


0.39 


0.22 


0.17 


0.22 


133 


0.41 


0.18 


0.23 


0.18 


161 


0.42 


0.17 


0.24 


0.17 


183 


0.42 


0.16 


0.26 


0.16 


189 


0.42 


0.16 


0.26 


0.16 



TABLE II: Branching ratios for the e'^e — > qq annihilations into various quark flavours as a function of centre- 
of-mass energy. 



The mass of resonances with F > 1 MeV has been distributed according to a relativistic Breit-Wigner 
function over the interval [mo — Ato,771o + Am], mo where mo is the central mass value and Am — 
min{2r,mo — mth}, T being the width and mth the threshold mass value for the opening of all allowed 
decay channels. The primary production rate of neutral mesons such as 77, 77', (j), lu and others, which 
are a superposition of ss and uu, dd states, has been suppressed with 7^ according to the fraction of ss 
content. For this purpose we use the mixing formulas quoted in the Review of Particles Properties [l^] 
with angles 9 = —10° and 6 = 39° for the (ry, 77') system and for (w, (f)) system, respectively, while for 
other nonets we used 9 = 28°. 

For each experiment, the most recent measurements have been considered. Multiple measurements 
from different experiments have been averaged according to the PDG method ^2^, with error rescaling 
in case of discrepancy, that is a /dof of the weighted average > 1. The overall calculated yields Ti are 
compared to the experimental measurements Ei, and the total overall y^, 

^^^Y.(T^-E^f/aj (24) 

i 

where at are the experimental errors, is minimized. The minimization is in fact carried out in two steps, 
in order to take into account the uncertainties on input parameters, such as hadron masses, widths and 
branching ratios, according to the following procedure [3|. First a with only experimental errors is 
minimised and preliminary best-fit model parameters are determined. Then, keeping the preliminarly 
fitted parameters fixed, the variations Ati^*''™ of the multiplicities corresponding to the variations of the 
Z*^ input parameter by one standard deviation are calculated. Such variations are considered as additional 
systematic uncertainties on the multiplicities and the following covariance matrix is formed. 
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(25) 
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to be added to the experimental covariance matrix C°^^. Finally a new is minimised with covariance 
matrix C'^p + C^'^^, from which the best-fit estimates of the parameters and their errors are obtained. 
Actually, more than 350 among the most relevant or poorly known input parameters have been varied. 
However, it should be mentioned that no variation of the branching fractions of heavy flavoured hadrons 
has been done. Therefore, for some specific species, the systematic error could have been underestimated. 



V. RESULTS 



A. Light flavoured hadrons 



We begin our analysis with the most extensive sample, the LEP data at 91.25 GeV. It comes from 
a compilation of results from the four different experimental groups (see references at the end of the 
Appendix), and it lists up to 30 different light flavoured species. However, for short-lived and hence 
broad resonances, the separation of resonance signal from background often becomes difficult, making 
the assessment of systematic errors problematic. Moreover, broad resonances yields are more sensitive 
to feeding from possibly unobserved heavier states or poorly known branching ratios. For this reason, 
we first consider, both for the conventional and for the Hawking-Unruh scenario, the analysis of the 
unproblematic ("golden") species of widths less than 10 MeV; this still leaves 15 different hadronic states 
to be analysed, listed in tables Hill and IIVI 
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FIG. 4: Comparison between measured and fitted multiplicities of long-lived hadronic species in e^e~ collisions 
at v^s = 91.25 GeV. Left: statistical hadronization model with one temperature. Right: Hawking-Unruh radiation 
model. 



As noted, the conventional statistical resonance gas approach is based on a universal temperature T, 
a strangeness suppression factor 7s, and a global volume V. The fit of the long-lived species is shown in 
detail in table IIIII and fig. [4] and the resulting fit parameters are 

r= 164.6 ±3.0 MeV; 7^ = 0.648 ± 0.026; = 40.2 ± 5.7 fm^ (26) 

with a xVdof = 39/12. The errors on the parameters are the fit errors rescaled by /dof- Such a 
method takes into account the additional uncertainty on the parameters if the fit leads to x^/dof > 1. 
This rescaling has been applied to all parameter errors quoted in this paper. 

Next, we perform the corresponding hadron-rcsonance gas analysis in the Hawking-Unruh formulation, 
introducing different temperatures determined by the string tension a and the strange quark mass mj. 
The results for long-lived species are shown in table IIVI and fig. 21 The resulting fit parameters here are 



a = 0.1683 ± 0.0048 GeV^; m, = 0.083 ± 0.004 GeV, 1/ = 40.3 ± 3.2 fm^; 



(27) 
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with a x^/dof = 22/12, somewhat better than that of the corresponding conventional fit. 

We now repeat both analyses using the entire 91.25 GeV data set, with the results shown in table XX 
and XXI of the Appendix. The resulting fit values (see tables IIIII and lIVp agree well within errors with 
those obtained from the "golden" data set at 91.25 GeV. As expected, because of the mentioned error 
sizes, the /dof for the full 91.25 set is considerably worse. 
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-4.98 


E+ 


0.0535 


± 


0.0052 


0.0429 


-2.0 


-19.9 


E° 


0.0389 


± 


0.0041 


0.0435 


1.1 


11.8 


E" 


0.0410 


± 


0.0037 


0.0391 


-0.51 


-4.58 




0.01319 


± 


0.00050 


0.01256 


-1.3 


-4.81 


Q. 


0.00062 


± 


0.00010 


0.00089 


2.7 


43.7 



TABLE III: Abundances of long-lived hadrons in e^e~ collisions at — 91.25 GeV, compared to a statistical 
hadronization model fit based on T and 7^. The third column shows the residual, defined as the difference 
between model and data divided by the error, while fourth column shows the differences between model and data 
in percent. References to the original experimental publications can be found in table IXXl in the Appendix. 





Experiment (E) 


Model (M) 


Residual 


IM - El/E [%] 




9.61 


± 


0.29 


9.73 


0.41 


1.25 


7r+ 


8.50 


± 


0.10 


8.32 


-1.7 


-2.06 


K+ 


1.127 


± 


0.026 


1.106 


-0.80 


-1.85 


K% 


1.0376 


± 


0.0096 


1.0656 


2.9 


2.69 


V 


1.059 


± 


0.086 


1.006 


-0.61 


-4.98 




1.024 


± 


0.059 


0.967 


-0.97 


-5.58 


P 


0.519 


± 


0.018 


0.559 


2.2 


7.78 




0.166 


± 


0.047 


0.093 


-1.6 


-43.8 


<!> 


0.0977 


± 


0.0058 


0.1057 


1.4 


8.11 


A 


0.1943 


± 


0.0038 


0.1892 


-1.3 


-2.63 


E+ 


0.0535 


± 


0.0052 


0.0438 


-1.9 


-18.2 


E" 


0.0389 


± 


0.0041 


0.0444 


1.4 


14.2 


E" 


0.0410 


± 


0.0037 


0.0401 


-0.25 


-2.22 




0.01319 


± 


0.00050 


0.01265 


-1.1 


-4.11 


Q. 


0.00062 


± 


0.00010 


0.00077 


1.5 


23.5 



TABLE IV: Abundances of long-lived hadrons in e"'"e~ collisions at -^s — 91.25 GeV, compared to a Hawking- 
Unruh fit in terms of string tension <j and strange quark mass ttIs. The third column shows the residual, defined 
as the difference between model and data divided by the error, while fourth column shows the differences between 
model and data in percent. References to the original experimental publications can be found in table IXXII in 
the Appendix 

Here a comment is in order. The simple formulae ^ and ([Tj), in both models, rely on some side 
assumptions (e.g. the special distributions for cluster charge fluctuations needed for the introduction of 
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the equivalent global cluster) that are not expected to be exactly fulfilled. Therefore, those formulae 
are to be taken as a zero-order approximation and not as a faithful representation of the real process. 
Deviations from the introduced assumption entail corrections to the formulae ([5]) and ([7]) which are 
nevertheless very difficult to estimate. The theoretical error involved in these formulae becomes important 
when the accuracy of measurements is comparable and, in this case, a bad is to be expected. This 
is probably the case at ^/s = 91.25 GeV, where the relative accuracy of measurements is of the order of 
few percent for many particles. In this case, the fit is a useful tool to determine the best parameters 
of the "simplified" theory but should be used very carefully as a measure of the fit quality. As has been 
mentioned, in order to take into account the uncertainty on parameters implied in fits with x^/do/ > 1, 
parameter errors have been rescaled by /dof if this is larger than 1, according to Particle Data Group 
procedure 

For all the remaining energies we have also carried out the corresponding analyses; the results are listed 
in tables |V] and IVII for the model parameters, while the comparison between measured and calculated 
multiplicities are shown in tables X to XXXI of the Appendix. 





T [MeV] 




7s 


xVdof 


14 


172.1 ± 5.2 


8.3 ± 1.0 


0.772 


± 


0.094 


0.9 / 3 


22 


178.7 ± 3.7 


8.70 ± 0.94 


0.76 


± 


0.10 


0.7 / 3 


29 


164.0 ± 5.4 


15.0 ± 2.4 


0.683 


± 


0.075 


33 / 13 


35 


163.3 ± 3.2 


15.0 ± 1.4 


0.730 


± 


0.045 


8.2 / 7 


43 


169 ± 10 


13.5 ± 3.2 


0.741 


± 


0.074 


2.9 / 3 


91 


161.9 ± 4.1 


25.8 ± 3.4 


0.638 


± 


0.039 


215 / 27 


91* 


164.6 ± 3.0 


23.3 ± 2.2 


0.648 


± 


0.026 


39 / 12 


133 


167.1 ± 7.5 


26.0 ± 4.6 


0.671 


± 


0.074 


0.1 / 2 


161 


153.4 ± 6.5 


37.2 ± 5.9 


0.72 


± 


0.12 


0.03 / 1 


183 


161 ± 13 


35 ± 11 


0.446 


± 


0.098 


5.0 / 2 


189 


159 ± 12 


36 ± 10 


0.54 


± 


0.11 


7.5 / 2 



TABLE V: Best fit parameters for the statistical hadronization model in e^e collisions. The golden sample fit 
is marked with a *. 





a [GcV^] 




[MeV] 


xVdof 


14 


0.185 


± 


0.015 


133 ± 24 


71 ± 19 


0.9 / 3 


22 


0.199 


± 


0.017 


140 ± 32 


77 ± 20 


0.7 / 3 


29 


0.1673 


± 


0.0096 


240 ± 36 


78 ± 15 


38 / 13 


35 


0.1675 


± 


0.0065 


237 ± 23 


74.6 ± 6.5 


8.8 / 7 


43 


0.178 


± 


0.020 


216 ± 48 


76 ± 16 


3.2 / 3 


91 


0.1625 


± 


0.0078 


406 ± 52 


82.3 ± 7.8 


217 / 27 


91* 


0.1683 


± 


0.0042 


368 ± 24 


83.2 ± 4.0 


23 / 12 


133 


0.175 


± 


0.015 


418 ± 69 


89 ± 16 


0.1 / 2 


161 


0.148 


± 


0.029 


590 ± 220 


74 ± 24 


0.03 / 1 


183 


0.165 


± 


0.026 


550 ± 160 


130 ± 28 


5.1 / 2 


189 


0.161 


± 


0.029 


560 ± 180 


110 ± 36 


7.7 / 2 



TABLE VI: Best fit parameters for the Hawking-Unruh model in e^e collisions. The golden sample fit is marked 
with a *. 



The parameters T, 7^ and cr, are shown in figs. [5] and [6] respectively. It can be seen that they are 
remarkably constant throughout the examined energy range from ^/s = 14 to 189 GeV. 

We now recall that the quantities we treated as fit parameters in the Hawking-Unruh analyses, the 
string tension and the strange quark mass, have in fact been determined in various other contexts and 
by different methods; they are quite well-known. The string tension a is obtained in studies of heavy 
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V s (GeV) 




V s (GeV) 



FIG. 5: Hadronization temperature T (left) and strangeness suppression factor 7s (right) from conventional 
statistical fits to hadron abundances in e^e~ annihilation, as function of the incident energy -/s. 
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FIG. 6: String tension (left) and strange quark mass (right) from Hawking-Unruh fits to hadron abundances in 
e^e~ annihilation, as function of the incident energy ^/s. The shaded bands give the overall average values as 
determined by other data. 



quarkonium spectroscopy as well as from Regge phenomenology. The canonical value [221 was 0.192 GeV^; 
more recent calculations range from 0.16 GeV^ [23] to 0.22 GeV^ f24j, giving an estimate of tr = 0.19±0.03. 
The best average value of the strange quark mass is presently listed as 20] nis ~ 0.095 ± 0.025 GeV. In 
both cases we have good agreement with our fit values, as seen in fig. [HI We can thus indeed conclude 
that the Hawking-Unruh approach provides a parameter-free description of thermal hadron abundances 
in e+e^ annihilation. The suppression of hadrons containing strange quarks is fully accounted for in 
terms of slight temperature changes due to the heavier strange quark mass. It is thus natural that this 
affects also non-strange hadrons dominantly made up of a strange quarks, such as the (j). 

To illustrate the effect, we list in table I VIII the different temperatures resulting from the different 
strange quark contents of the observed hadrons, using the fit parameters from the golden 91.25 data set. 
To check what such temperature differences can lead to, we compare the rate of direct production in 



the conventional to that of the Hawking-Unruh scenario. This direct rate is given by 
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(^^>0 = 3^^K2(m/T)7§ (28) 

in the conventional scenario; using the values ([26|) together with the production value listed in table |Vl 
we obtain (n)^ ~ 0.078. Note that 7^ ~ 0.42 reduces the equilibrium value by more than a factor of two. 
The Hawking-Unruh scheme has with T{ss) — 0.148 GeV a lower temperature for a meson containing a 
strange quark-antiquark pair than that governing light quark mesons. With 

(">0 = 3S£^K2(m/r(s.)) (29) 

and the corresponding production volume of table IVTl this leads to (ri)^ ~ 0.077 and hence practically 
the same value, however, without invoking the parameter 7^. We note that these results should not be 
compared directly to the </> production measured in e'^e^ annihilation, which contains (at 91.25 GeV) a 
further 30 % due to feed-down contributions from charmed and bottomed hadron decay. 



T [GeV] 




T(00) 


0.164 


T{Os) 


0.156 


Tiss) 


0.148 


T(OOO) 


0.164 


T(OOs) 


0.158 


r(Oss) 


0.153 


T{sss) 


0.148 



TABLE VII: Hadronization temperatures for hadrons of different strangeness content, for — 0.083 Gev and 
a = 0.169 GeV^ 



B. Heavy flavoured hadrons 

As has been mentioned in the previous section, the calculation of heavy flavoured hadron yields in 
e'^e~ —^ cc and e~^e~ — > bb events is necessary to determine the final light flavoured hadron multiplicities. 
The heavy flavoured hadron relative abundances in such events are determined according to formula ([7])- 
For the conventional formulation of the statistical model, one has the same temperature as for the light- 
flavoured hadron species, while for the Hawking- Unruh radiation model, we do not have a deflnite 
prescription and we chose to use the same temperatures as for the light flavoured hadrons, which are 
dependent on the quark content. 

Once the model parameters have been fitted by using light-flavoured hadronic multiplicities, they can 
be used to predict relative yields of heavy flavoured species in e+e~ — > cc and e+e" bb annihilations 
and compare them to measured ones. This is a powerful, parameter-free, independent test of the con- 
ventional statistical model and a necessary consistency check for the Hawking-Unruh radiation model. 
Also, comparing theoretical values to measured relative abundances of heavy flavoured hadrons in speciflc 
annihilations channels (e.g. charmed hadrons in e+e~ cc) we achieve a more effective test because, 
e.g., the contribution of weak b ^ c decays is excluded. 

The relative yields of several heavy flavoured hadronic species have been measured in e+e" collisions at 
91.25 GeV by the four LEP experiments. We show a comparison between model and weighted averaged 
experimental values in tables IVIIII and IIXI 

For the statistical model (table [VnTI, the theoretical values have been estimated by using the param- 
eters in eq. (j26p . The agreement between model and experiment is strikingly good, with few peculiar 
deviations in heavier states. In all those cases we observe an underestimation of measured values which 
may partly be explained by the absence, in our input spectrum, of unknown heavier heavy flavoured 
resonances feeding these states. It is quite remarkable that the model is able to reproduce the largely 
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Particle 




Experiment (E) 


Model (M) 


Residual 


[M-E)/E [%] 


D« 


[28J 


0.559 ± 0.022 


0.5406 


-0.83 


-3.2 


D+ 


[28] 


0.238 ± 0.024 


0.2235 


-0.60 


-6.1 


D*+ 


[28, 29, 30] 


0.2377± 0.0098 


0.2279 


-1.00 


-4.1 


D*° 


m 


0.218 ± 0.071 


0.2311 


0.18 


6.0 


Dl 


[32, 33] 


0.0173 ± 0.0039 


0.01830 


0.26 


5.8 


Df 


[32, 33] 


0.0484 ± 0.0080 


0.02489 


-2.94 


-48.6 




[28] 


0.116 ± 0.036 


0.1162 


0.006 


0.19 


Dl 


[28] 


0.069 ± 0.026 


0.0674 


-0.06 


-2.4 


Dsi 


[33, 34] 


0.0106 ± 0.0025 


0.00575 


-1.94 


-45.7 


DI2 


[3^ 


0.0140 ± 0.0062 


0.00778 


-1.00 


-44.5 


Ac 


[28] 


0.079 ± 0.022 


0.0966 


0.80 


22.2 


{B° + B+)/2 


[35] 


0.399 ± 0.011 


0.3971 


-0.18 


-0.49 


Bs 


[35] 


0.098 ± 0.012 


0.1084 


0.87 


10.6 


ByB{nda) 


[36, 37, 38, 39] 


0.749 ± 0.040 


0.6943 


-1.37 


-7.3 


B" X BR{B{*)ti) 


[40,41,42] 


0.180 ± 0.025 


0.1319 


-1.92 


-26.7 


(B^ + Bi) X BR{B{*)tv) 


[4^ 


0.090 ± 0.018 


0.0800 


-0.57 


-11.4 


BI2 X BR{BK) 


[4^ 


0.0093 ± 0.0024 


0.00631 


-1.24 


-32.1 


b-baryon 


[35] 


0.103 ± 0.018 


0.09751 


-0.30 


-5.3 


"6 


[35] 


0.011 ± 0.006 


0.00944 


-0.26 


-14.2 



TABLE VIII: Abundances of charmed hadrons in e^e cc annihilations and bottomed hadrons in e^e 
bb annihilations at ^/s = 91.25 GeV, compared to the prediction of the statistical model. 



different V/P ratios in the charm {D* /D) and bottomed {B* / B) sector - a long-standingissue in string 
models - without any additional parameter. This result confirms previous early findings TH. . 

For the Hawking- Unruh radiation model (table HXj) . the agreement between theoretical multiplicities, 
calculated with the parameters in eq. (jTZ]) and experiment is generally good, although not as good as for 
the conventional scheme. In fact, there are some specific discrepancies, especially in the beauty sector. In 
particular, we underestimate the relative yield of Bs mesons, which is an effect of the lower temperature 
for open strange particles combined with the high mass of these particles, compared to light-flavoured 
species. In general, the way heavy flavoured hadrons are to be calculated in this model is still an open 
issue. 



VI. CONCLUSIONS 

We have shown that in accord with previous studies [H, Q, Q , the thermal hadron abundances observed 
in e''"e~ collisions over a wide range of energies can indeed be accounted for in an ideal resonance gas 
scenario, based on a universal temperature T ~ 165 MeV and a strangeness suppression factor 7s ~ 0.7. 
The latter is the ad hoc price paid in order to account for the deviation from full chemical equilibrium 
observed in the data. Remarkably, also the relative abundances of heavy flavoured species are in very 
good agreement with the statistical-thermal ansatz. 

The Hawking-Unruh scenario, on the other hand, provides an intrinsic deviation from full equilibrium 
through the dependence of the radiation temperature on the mass of the emitted quark. Given the value 
of this mass and the string tension specifying the fleld strength at the conflnement horizon, we then have 
a parameter-free prediction of the relative hadron abundances. We have seen here that these predictions 
agree well with the data at all energies, with the caveat that the relative multiplicities of heavy flavoured 
hadron species are not completely understood in this picture and are in slightly worse agreement with 
respect to the conventional statistical model. In a subsequent paper, we shall extend the Hawking-Unruh 
description to high energy heavy ion collisions, where it becomes significantly modified. 
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Particle 




Experiment (E) 


Model (M) 


Residual 


{M-E)/E [%] 


D« 


[28J 


0.559 ± 0.022 


0.5635 


0.20 


0.81 


D+ 


[28] 


0.238 ± 0.024 


0.2332 


-0.20 


-2.0 


D*+ 


[28, 29, 30] 


0.2377± 0.0098 


0.2373 


-0.04 


-0.18 


D*° 


m 


0.218 ± 0.071 


0.2407 


0.32 


10.4 


Dl 


[32, 33] 


0.0173 ± 0.0039 


0.01897 


0.43 


9.6 


Df 


[32, 33] 


0.0484 ± 0.0080 


0.02577 


-2.82 


-46.8 




[28] 


0.116 ± 0.036 


0.08460 


-0.87 


-27.1 


Dl 


[28] 


0.069 ± 0.026 


0.04793 


-0.81 


-30.5 


Dsi 


[33, 34] 


0.0106 ± 0.0025 


0.00356 


-2.82 


-66.5 


DI2 


[3^ 


0.0140 ± 0.0062 


0.00479 


-1.49 


-66.1 


Ac 


[28] 


0.079 ± 0.022 


0.09922 


0.92 


25.6 


{B° + B+)/2 


[35] 


0.399 ± 0.011 


0.4390 


3.64 


10.0 


Bs 


[35] 


0.098 ± 0.012 


0.0276 


-5.87 


-71.9 


ByB{nda) 


[36, 37, 38, 39] 


0.749 ± 0.040 


0.6978 


-1.28 


-6.8 


B" X BR{B{*)ti) 


[40, 41, 42] 


0.180 ± 0.025 


0.1479 


-1.28 


-17.8 


(B^ + Bi) X BR{B{*)tv) 


[4^ 


0.090 ± 0.018 


0.0894 


-0.04 


-0.72 


BI2 X BR{BK) 


[4^ 


0.0093 ± 0.0024 


0.00136 


-3.31 


-85.3 


b-baryon 


[35] 


0.103 ± 0.018 


0.0944 


-0.48 


-8.4 


"6 


[35] 


0.011 ± 0.006 


0.00415 


-1.14 


-62.2 



TABLE IX: Abundances of charmed hadrons in e'''e —^ cc annihilations and bottomed hadrons in e'^e 
bb annihilations at y/s = 91.25 GeV, compared to the prediction of the Hawking-Unruh radiation model. 



In closing, we comment on the degree of agreement between experiment and theory in our description. 
Hawking-Unruh radiation is thermal in leading order [2^ . with higher order interaction terms. Similarly, 
one expects corrections to the simplest statistical hadronization model formulae, see discussion in Sect. 4. 
When accuracy of measurements is good enough, such higher-order effects must be taken into account 
and the fit quality to the simplest formulae unavoidably degrades. It is an interesting question to see if 
nuclear collisions, with a higher degree of averaging, lead to smaller deviations with measurements of the 
same accuracy. 

Note added in proof 

After completion of this work, another analysis of hadron production in e+e^ annihilation has appeared, 
reaching very different conclusions [27| . However, in contrast to our work, it is assumed in that analysis 
that the conservation of charm and bottom can be neglected. In e+e^ annihilation, more than 40% of 
all events contain a primary charm or bottom quark- antiquark pair (see our table |ll|, hence two heavy- 
flavoured hadrons. Neglecting the corresponding conservation conditions and the decay contributions 
of these heavy-flavoured hadrons into light-flavoured hadrons necessarily produces serious disagreement 
with the data, especially for strange particles. 
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Appendix 

In the following Tables |X] - IXXIXl the experimental data and the statistical hadronization model 
predictions for various hadron multiplicities are compared. The first column shows the experimental 
value while in the second column, the statistical hadronization model prediction is quoted. In the fourth 
column one can find the relative deviation of the model from the experiment in percentage while the 
third column shows the residuals defined as 

Mth _ jyex 

Residuali = ^ (30) 

in which Nf^ and Nf^ are the theoretical and experimental multiplicitites and ai is the (experimental) 
standard error of a particle species i. For each of the energies, we show 2 different tables one after the 
other so that always the first table shows the conventional statistical hadronization model fit results while 
the following one shows the same information in the Hawking-Unruh approach. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[43] 


4.69 ± 0.20 


4.65 


-0.18 


-0.752 


7r+ 




3.60 ± 0.30 


3.79 


0.62 


5.18 


K+ 


m 


0.600 ± 0.070 


0.589 


-0.15 


-1.81 


K% 


[44, 45] 


0.563 ± 0.045 


0.556 


-0.15 


-1.17 


P 


m 


0.210 ± 0.030 


0.199 


-0.37 


-5.36 


A 


m 


0.065 ± 0.020 


0.077 


0.58 


17.8 



TABLE X: Hadron multiplicities in e^e~ collisions at 14 GeV, compared to the outcomes of the fit based on 
the statistical hadronization model with T and 75. The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[43] 


4.69 ± 0.20 


4.66 


-0.17 


-0.712 


Tr+ 


m 


3.60 ± 0.30 


3.79 


0.62 


5.16 


K+ 


m 


0.600 ± 0.070 


0.589 


-0.16 


-1.85 


K% 


[44,45] 


0.563 ± 0.045 


0.556 


-0.14 


-1.14 


P 


[44] 


0.210 ± 0.030 


0.198 


-0.39 


-5.52 


A 


[44] 


0.065 ± 0.020 


0.077 


0.58 


17.8 



TABLE XI; Hadron multiplicities in e'''e~ collisions at 14 GeV, compared to the outcomes of the fit based on 
Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model in 
units of standard error and in percentages, respectively. 
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Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[43] 


5.50 ± 0.40 


5.49 


-0.033 


-0.238 


n+ 


m] 


4.40 ± 0.50 


4.54 


0.28 


3.16 


K+ 


m 


0.75 ± 0.10 


0.69 


-0.65 


-8.64 


K% 


[45,46] 


0.638 ± 0.057 


0.651 


0.22 


1.96 


P 


[44] 


0.310 ± 0.030 


0.305 


-0.16 


-1.51 


A 


[44] 


0.110 ± 0.025 


0.117 


0.29 


6.70 



TABLE XII: Hadron multiplicities in e"'"e~ collisions at 22 GeV, compared to the outcomes of the fit based on 
the statistical hadronization model with T and js- The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[43] 


5.50 ± 0.40 


5.49 


-0.020 


-0.148 


n+ 


m 


4.40 ± 0.50 


4.54 


0.28 


3.18 


K+ 


m 


0.75 ± 0.10 


0.68 


-0.65 


-8.71 


K% 


[45,46] 


0.638 ± 0.057 


0.651 


0.22 


1.94 


P 


m 


0.310 ± 0.030 


0.305 


-0.17 


-1.65 


A 


[44] 


0.110 ± 0.025 


0.118 


0.31 


7.02 



TABLE XIII: Hadron multiplicities in e'''e~ collisions at 22 GeV, compared to the outcomes of the fit based on 
Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model in 
units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[47] 


5.30 


± 


0.70 


6.48 


1.7 


22.2 


7r+ 


[48] 


5.35 


± 


0.25 


5.42 


0.26 


1.23 


K+ 


[48] 


0.700 


± 


0.050 


0.747 


0.93 


6.66 


K% 


[49, 50, 5^ 


0.691 


± 


0.029 


0.712 


0.73 


3.05 


V 


[52, 53] 


0.584 


± 


0.075 


0.654 


0.92 


11.8 


p' 


M 


0.900 


± 


0.050 


0.745 


-3.1 


-17.2 


K*+ 


[55] 


0.310 


± 


0.030 


0.237 


-2.4 


-23.7 


K*° 


[50, 5^ 


0.281 


± 


0.022 


0.232 


-2.2 


-17.3 


P 


[48] 


0.300 


± 


0.050 


0.300 


0.0038 


0.0627 


v' 


[53] 


0.26 


± 


0.10 


0.07 


-1.8 


-71.6 






0.084 


± 


0.022 


0.092 


0.35 


9.15 


A 


[57, 58, 59] 


0.0983 


± 


0.0060 


0.1016 


0.56 


3.43 




[60, ei 


0.0083 


± 


0.0020 


0.0070 


-0.64 


-15.4 


E*+ 


m 


0.0083 


± 


0.0024 


0.0111 


1.2 


34.2 


k;+ 


[55] 


0.045 


± 


0.022 


0.016 


-1.3 


-64.4 


n 


m 


0.0070 


± 


0.0036 


0.0005 


-1.8 


-93.4 



TABLE XIV: Hadron multiplicities in e^e~ collisions at 29 GeV, compared to the outcomes of the fit based on 
the statistical hadronization model with T and 75. The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 
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Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[47] 


5.30 


± 


0.70 


6.37 


1.5 


20.2 


7r+ 


[48] 


5.35 


± 


0.25 


5.31 


-0.15 


-0.715 


K+ 


[48] 


0.700 


± 


0.050 


0.760 


1.2 


8.55 


K% 


[49, 50, 5^ 


0.691 


± 


0.029 


0.725 


1.2 


4.96 


V 


[52, 53] 


0.584 


± 


0.075 


0.643 


0.78 


10.0 


P° 


[5^ 


0.900 


± 


0.050 


0.727 


-3.5 


-19.2 


K*+ 


[55] 


0.310 


± 


0.030 


0.231 


-2.6 


-25.6 


K*° 


[50, 5^ 


0.281 


± 


0.022 


0.227 


-2.4 


-19.3 


P 


[48] 


0.300 


± 


0.050 


0.292 


-0.17 


-2.82 


v' 


[53] 


0.26 


± 


0.10 


0.07 


-1.9 


-75.0 




[56] 


0.084 


± 


0.022 


0.084 


-0.022 


-0.593 


A 


[57, 58, 59] 


0.0983 


± 


0.0060 


0.1024 


0.69 


4.24 


H- 


[60, 61J 


0.0083 


± 


0.0020 


0.0068 


-0.72 


-17.5 


E*+ 


[6^ 


0.0083 


± 


0.0024 


0.0111 


1.2 


34.8 


K*+ 


[55] 


0.045 


± 


0.022 


0.014 


-1.4 


-69.9 


n 


[62] 


0.0070 


± 


0.0036 


0.0004 


-1.8 


-94.6 



TABLE XV: Hadron multiplicities in e"''e~ collisions at 29 GeV,compared to the outcomes of the fit based on 
Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model in 
units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[63, 64, 65] 


6.31 ± 0.35 


6.48 


0.49 


2.73 


7r+ 


[66] 


5.45 ± 0.25 


5.42 


-0.14 


-0.621 


K+ 


[66] 


0.88 ± 0.10 


0.78 


-0.98 


-11.2 


K% 


[45, 46, 67] 


0.740 ± 0.017 


0.746 


0.33 


0.759 


V 


[64, 65] 


0.636 ± 0.080 


0.661 


0.32 


4.06 


P° 


[44, 68] 


0.756 ± 0.077 


0.739 


-0.23 


-2.30 


K*+ 


[46, 67, 68] 


0.361 ± 0.046 


0.248 


-2.4 


-31.2 


P 


[66, 69] 


0.302 ± 0.033 


0.300 


-0.078 


-0.838 


A 


[67, 70] 


0.108 ± 0.010 


0.108 


-0.042 


-0.391 




[70] 


0.0060 ± 0.0021 


0.0079 


0.90 


31.5 



TABLE XVL Hadron multiplicities in e^e~ collisions at 35 GeV, compared to the outcomes of the fit based on 
the statistical hadronization model with T and 75. The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 
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Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[63, 64, 65) 


6.31 ± 0.35 


6.46 


0.41 


2.31 


7V + 


[66j 


5.45 ± 0.25 


5.39 


-0.25 


-1.13 


K+ 


[66] 


0.88 ± 0.10 


0.78 


-0.97 


-11.0 


K% 


[45, 46, 67) 


0.740 ± 0.017 


0.748 


0.47 


1.08 


V 


[64,65] 


0.636 ± 0.080 


0.657 


0.26 


3.33 


P° 


[44, 68] 


0.756 ± 0.077 


0.735 


-0.28 


-2.82 


K*+ 


[46, 67, 68] 


0.361 ± 0.046 


0.239 


-2.6 


-33.6 


P 


[66, 69] 


0.302 ± 0.033 


0.301 


-0.042 


-0.450 


A 


[67, 70] 


0.108 ± 0.010 


0.108 


0.028 


0.261 




[70] 


0.0060 ± 0.0021 


0.0075 


0.70 


24.4 



TABLE XVII: Hadron multiplicities in e'''e~ collisions at 35 GeV,compared to the outcomes of the fit based on 
Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model in 
units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[64, 66] 


6.66 ± 0.65 


6.63 


-0.055 


-0.541 


n+ 


[66] 


5.55 ± 0.25 


5.55 


-0.0021 


-0.00955 


K+ 


m 


0.96 ± 0.15 


0.81 


-1.0 


-15.9 


K's 




0.760 ± 0.035 


0.775 


0.43 


2.01 


K*+ 


[46] 


0.385 ± 0.094 


0.264 


-1.3 


-31.3 


A 


[70] 


0.128 ± 0.024 


0.131 


0.12 


2.23 



TABLE XVIII: Hadron multiplicities in e^e~ collisions at 43 GeV, compared to the outcomes of the fit based on 
the statistical hadronization model with T and 7s. The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[64, 66] 


6.66 ± 0.65 


6.62 


-0.065 


-0.639 


n+ 


[66] 


5.55 ± 0.25 


5.54 


-0.042 


-0.188 


K+ 


[71] 


0.96 ± 0.15 


0.81 


-1.0 


-15.8 


Kl 


[46] 


0.760 ± 0.035 


0.777 


0.47 


2.18 


K*+ 


[46] 


0.385 ± 0.094 


0.255 


-1.4 


-33.6 


A 


[70] 


0.128 ± 0.024 


0.131 


0.13 


2.52 



TABLE XIX: Hadron multiplicities in e^e~ collisions at 43 GeV, compared to the outcomes of the fit based on 
Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model in 
units of standard error and in percentages, respectively. 
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Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[72, 73, 74, 75] 


9.61 


± 


0.29 


10.20 


2.0 


6.18 


7r+ 


[76, 77, 78, 79] 


8.50 


± 


0.10 


8.76 


2.6 


3.13 


K+ 


[76, 77, 78, 79] 


1.127 


± 


0.026 


1.091 


-1.4 


-3.20 


Kl 


[72, 75, 79, 80, 83) 


1.0376 


± 


0.0096 


1.0507 


1.4 


1.26 


V 


[72, 75, 8^ 


1.059 


± 


0.086 


1.041 


-0.21 


-1.70 


P° 


[85, 86] 


1.40 


± 


0.13 


1.19 


-1.6 


-15.0 


P+ 


[sa 


1.20 


± 


0.22 


1.14 


-0.26 


-4.66 




[84,87, 88] 


1.024 


± 


0.059 


1.014 


-0.17 


-0.997 


K*+ 


[80, 82, 89] 


0.357 


± 


0.022 


0.353 


-0.16 


-1.02 




[79, 82, 90, 91J 


0.370 


± 


0.013 


0.346 


-1.9 


-6.33 


P 


[76, 77, 78, 79] 


0.519 


± 


0.018 


0.564 


2.5 


8.76 


v' 


[87, 88] 


0.166 


± 


0.047 


0.106 


-1.3 


-36.1 


/o 


[85, 86, 92] 


0.1555 


± 


0.0085 


0.0779 


-9.1 


-49.9 


a-o 


[82 


0.135 


± 


0.054 


0.084 


-0.95 


-37.8 


4> 


[79, 82, 91, 92] 


0.0977 


± 


0.0058 


0.1150 


3.0 


17.7 


A 


[72, 79, 83, 93, 94] 


0.1943 


± 


0.0038 


0.1779 


-4.3 


-8.42 


E+ 


[95, 96, 97] 


0.0535 


± 


0.0052 


0.0415 


-2.3 


-22.4 


T° 


[82, 96, 97, 98] 


0.0389 


± 


0.0041 


0.0421 


0.77 


8.11 


E" 


[96, 99] 


0.0410 


± 


0.0037 


0.0378 


-0.85 


-7.65 


A++ 


[100, io^ 


0.044 


± 


0.017 


0.090 


2.7 


105. 


/2 


[85, 86, 92] 


0.188 


± 


0.020 


0.122 


-3.4 


-35.1 


/l 


[102] 


0.165 


± 


0.051 


0.064 


-2.0 


-61.5 




[82,94, 95] 


0.01319 


± 


0.00050 


0.01187 


-2.6 


-10.0 


E*+ 


[82,94, 95] 


0.0118 


± 


0.0011 


0.0201 


7.5 


70.3 




[102] 


0.056 


± 


0.012 


0.010 


-3.9 


-82.6 




[85] 


0.036 


± 


0.012 


0.026 


-0.91 


-29.2 


A(1520) 


[94, 99] 


0.0112 


± 


0.0014 


0.0109 


-0.22 


-2.73 


/2 


[85] 


0.0120 


± 


0.0058 


0.0100 


-0.34 


-16.6 


=■•(' 


[82,94, 95] 


0.00289 


± 


0.00050 


0.00417 


2.6 


44.4 




[82,94, 98] 


0.00062 


± 


0.00010 


0.00081 


1.9 


31.0 



TABLE XX: Hadron multiplicities in e"^e~ collisions at 91 GeV, compared to the outcomes of the fit based on 
the statistical hadronization model with T and 7s. The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 
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Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[72, 73, 74, 75] 


9.61 


± 


0.29 


10.01 


1.4 


4.18 


7r+ 


[76, 77, 78, 79] 


8.50 


± 


0.10 


8.57 


0.77 


0.918 


K+ 


[76, 77, 78, 79] 


1.127 


± 


0.026 


1.131 


0.16 


0.373 


Kl 


[72, 75, 79, 80, 83) 


1.0376 


± 


0.0096 


1.0901 


5.5 


5.05 


V 


[72, 75, 8^ 


1.059 


± 


0.086 


1.026 


-0.38 


-3.12 


P° 


[85, 86] 


1.40 


± 


0.13 


1.15 


-1.9 


-17.4 


P+ 


[sa 


1.20 


± 


0.22 


1.11 


-0.42 


-7.53 




[84,87, 88] 


1.024 


± 


0.059 


0.982 


-0.71 


-4.08 


K*+ 


[80, 82, 89] 


0.357 


± 


0.022 


0.345 


-0.54 


-3.36 




[79, 82, 90, 91J 


0.370 


± 


0.013 


0.338 


-2.5 


-8.62 


P 


[76, 77, 78, 79] 


0.519 


± 


0.018 


0.548 


1.6 


5.67 


v' 


[87, 88] 


0.166 


± 


0.047 


0.093 


-1.6 


-43.8 


/o 


[85, 86, 92] 


0.1555 


± 


0.0085 


0.0751 


-9.5 


-51.7 


a-o 


[82 


0.135 


± 


0.054 


0.081 


-1.0 


-40.0 


4> 


[79, 82, 91, 92] 


0.0977 


± 


0.0058 


0.1048 


1.2 


7.19 


A 


[72, 79, 83, 93, 94] 


0.1943 


± 


0.0038 


0.1826 


-3.1 


-6.04 


E+ 


[95, 96, 97] 


0.0535 


± 


0.0052 


0.0424 


-2.1 


-20.7 




[82, 96, 97, 98] 


0.0389 


± 


0.0041 


0.0430 


1.00 


10.5 


E" 


[96, 99] 


0.0410 


± 


0.0037 


0.0388 


-0.59 


-5.31 


A++ 


[100, io^ 


0.044 


± 


0.017 


0.086 


2.5 


95.0 


/2 


[85, 86, 92] 


0.188 


± 


0.020 


0.115 


-3.7 


-38.9 


/l 


[102] 


0.165 


± 


0.051 


0.061 


-2.0 


-63.2 




[82,94, 95] 


0.01319 


± 


0.00050 


0.01204 


-2.3 


-8.72 


E*+ 


[82,94, 95] 


0.0118 


± 


0.0011 


0.0204 


7.8 


72.8 




[102] 


0.056 


± 


0.012 


0.007 


-4.1 


-87.3 




[85] 


0.036 


± 


0.012 


0.021 


-1.3 


-41.7 


A(1520) 


[94, 99] 


0.0112 


± 


0.0014 


0.0106 


-0.45 


-5.55 


/2 


[85] 


0.0120 


± 


0.0058 


0.0068 


-0.89 


-43.1 


=■•(' 


[82,94, 95] 


0.00289 


± 


0.00050 


0.00423 


2.7 


46.4 




[82,94, 98] 


0.00062 


± 


0.00010 


0.00071 


0.85 


13.8 



TABLE XXI: Hadron multiplicities in e^e~ collisions at 91 GeV, compared to the outcomes of the fit based on 
Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model in 
units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M 


- E)/E (%) 


7r+ 


[103] 


9.92 ± 0.26 


9.94 


0.063 




0.167 


K+ 




1.30 ± 0.15 


1.29 


-0.064 




-0.714 


K% 


[103] 


1.25 ± 0.12 


1.25 


-0.070 




-0.663 


P 


[103] 


0.78 ± 0.13 


0.75 


-0.19 




-3.26 


A 


[103] 


0.250 ± 0.038 


0.256 


0.15 




2.27 



TABLE XXIL Hadron multiplicities in e^e~ collisions at 133 GeV, compared to the outcomes of the fit based on 
the statistical hadronization model with T and 75. The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 
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Experiment (E) 


Model (M) 


Residual 


(M 


- E)/E (%) 


n+ 


\m] 


9.92 ± 0.26 


9.94 


0.059 




0.157 


K+ 


[103] 


1.30 ± 0.15 


1.29 


-0.060 




-0.669 


Kl 


[103] 


1.25 ± 0.12 


1.25 


-0.056 




-0.533 


P 


[103] 


0.78 ± 0.13 


0.76 


-0.17 




-2.92 


A 


[103] 


0.250 ± 0.038 


0.255 


0.13 




1.92 



TABLE XXIII: Hadron multiplicities in e^e~ collisions at 133 GeV, compared to the outcomes of the fit based 
on Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model 
in units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 


n+ 


[103] 


10.38 ± 0.38 


10.37 


-0.00096 


-0.00355 


K+ 


[103] 


1.44 ± 0.30 


1.39 


-0.15 


-3.12 


K°s 


im 


1.32 ± 0.18 


1.34 


0.094 


1.31 


P 


[103] 


0.60 ± 0.24 


0.60 


-0.0010 


-0.0414 



TABLE XXIV: Hadron multiplicities in e^e~ collisions at 161 GeV, compared to the outcomes of the fit based 
on the statistical hadronization model with T and 75. The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 


7r+ 


[103] 


10.38 ± 0.38 


10.37 


-0.0012 


-0.00431 


K+ 


[103] 


1.44 ± 0.30 


1.39 


-0.15 


-3.18 


Kl 


[103] 


1.32 ± 0.18 


1.34 


0.097 


1.34 


P 


[103] 


0.60 ± 0.24 


0.60 


-0.0011 


-0.0436 



TABLE XXV: Hadron multiplicities in e+e^ collisions at 161 GeV, compared to the outcomes of the fit based on 
Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model in 
units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 




[103] 


10.89 ± 0.29 


10.87 


-0.081 


-0.216 


K+ 


\m 


1.42 ± 0.20 


1.03 


-2.0 


-27.2 


Kl 


[103] 


0.905 ± 0.086 


0.995 


1.0 


9.97 


P 


[103] 


0.66 ± 0.19 


0.71 


0.25 


7.13 


A 


[103] 


0.165 ± 0.025 


0.161 


-0.15 


-2.35 



TABLE XXVI: Hadron multiplicities in e"^e~ collisions at 183 GeV, compared to the outcomes of the fit based 
on the statistical hadronization model with T and 75. The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 


7r+ 


[103] 


10.89 ± 0.29 


10.86 


-0.12 


-0.326 


K+ 


[103] 


1.42 ± 0.20 


1.03 


-1.9 


-27.0 


K% 


[103] 


0.905 ± 0.086 


0.999 


1.1 


10.4 


P 


[103] 


0.66 ± 0.19 


0.73 


0.36 


10.5 


A 


[103] 


0.165 ± 0.025 


0.160 


-0.20 


-3.08 



TABLE XXVII: Hadron multiplicities in e'''e~ collisions at 183 GeV, compared to the outcomes of the fit based 
on Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model 
in units of standard error and in percentages, respectively. 



23 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 


n+ 




11.10 ± 0.26 


11.06 


-0.15 


-0.344 


K+ 


rio3] 


1.57 ± 0.16 


1.21 


-2.3 


-23.2 


Kl 


fios] 


1.060 ± 0.078 


1.169 


1.4 


10.3 


P 


fl03] 


0.59 ± 0.22 


0.72 


0.54 


20.3 


A 


fl03] 


0.200 ± 0.021 


0.196 


-0.20 


-2.16 



TABLE XXVIII: Hadron multiplicities in e"^e~ collisions at 189 GeV, compared to the outcomes of the fit based 
on the statistical hadronization model with T and 7s. The third and fourth columns show the differences between 
data and model in units of standard error and in percentages, respectively. 







Experiment (E) 


Model (M) 


Residual 


(M - E)/E (%) 


n+ 


[103] 


11.10 ± 0.26 


11.05 


-0.18 


-0.422 


K+ 


[103] 


1.57 ± 0.16 


1.21 


-2.3 


-23.2 


Kl 


[103] 


1.060 ± 0.078 


1.171 


1.4 


10.5 


P 


[103] 


0.59 ± 0.22 


0.74 


0.64 


23.7 


A 


fl03] 


0.200 ± 0.021 


0.195 


-0.22 


-2.31 



TABLE XXIX: Hadron multiplicities in e^e~ collisions at 189 GeV,compared to the outcomes of the fit based 
on Hawking-Unruh radiation model. The third and fourth columns show the differences between data and model 
in units of standard error and in percentages, respectively. 
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